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ABSTRACT 

We present the first large-scale effort of creating composite spectra of high-redshift type la super- 
novae (SNe la) and comparing them to low-redshift counterparts. Through the ESSENCE project, 
we have obtained 107 spectra of 88 high-redshift SNe la with excellent light-curve information. In 
addition, we have obtained 397 spectra of low-redshift SNe through a multiple-decade effort at Lick 
and Keck Observatories, and we have used 45 ultraviolet spectra obtained by HST/IUE. The low- 
redshift spectra act as a control sample when comparing to the ESSENCE spectra. In all instances, 
the ESSENCE and Lick composite spectra appear very similar. The addition of galaxy light to the 
Lick composite spectra allows a nearly perfect match of the overall spectral-energy distribution with 
the ESSENCE composite spectra, indicating that the high-redshift SNe are more contaminated with 
host-galaxy light than their low-redshift counterparts. This is caused by observing objects at all red- 
shifts with similar slit widths, which corresponds to different projected distances. After correcting 
for the galaxy-light contamination, subtle differences in the spectra remain. We have estimated the 
systematic errors when using current spectral templates for K-corrections to be ~0.02 mag. The 
variance in the composite spectra give an estimate of the intrinsic variance in low-redshift maximum- 
light SN spectra of ^3% in the optical and growing toward the ultraviolet. The difference between 
the maximum-light low and high-redshift spectra constrain SN evolution between our samples to be 
< 10% in the rest-frame optical. 

Subject headings: supernovae: general — cosmology: observations — cosmology — distance scale 
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Type la supernovae (SNe la) are the most precise 
known distance indicators at cosmological redshifts. The 
meticulous measurement of several hundred SNe la at 
both low and high redshifts has shown that the expan- 
sion of the Universe is currently acce l erating (iRiess et al.l 
1998', '2007*; 'Perlmutter.eLaI| ll999l; 'Astier e t al.l l200d 
Wood-Vascy et al. 2007); see Eirppenko (2G0£l) for a re- 
cent review. The underlying assumption behind that 
work is that high-redshift SNe la have the same peak lu- 
minosity as low-redshift SNe la (after corrections based 
on light-curve shape; e.g.. iPhillipsll 1 9931 ) . The luminosity 
of a given SN and its light-curve shape are determined 
by initial conditions of the white dwarf progenitor star 
(e.g., mass at explosion, C/0 abundance, and metallic- 
ity), and the properties of the explosion (e.g., deflagra- 
tion/detonation transition, the amount of unburned ma- 
terial, and the density at the ignition point). The pro- 
genitor properties are set by the initial conditions at the 
formation of the progenitor system, presumably having 
properties similar to the global galactic properties at that 
time. Since low-redshift SN la progenitor systems likely 
form, on average, in significantly different environments 
than high-redshift SN la progenitors, one may assume 
that some amount of evolution is inevitable (for a dis- 
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cus sion of difFerent c auses and effects of SN evolution, 
see lLeibundgutl[200lh . 

Theoretical studies of SN la evolution have focused on 
the composition, particularly metallicity, of the progen- 
itor system as the primary potential difference between 
the two samples. There have been two major studies 
with c onflicting results. For their study, iHoflich et al.l 
()1998l ) changed the progenitor metallicity and modeled 
the explosion, incl uding a full nuclear-reaction network. 
iLentz et al.l ([2000) chan ged th e metallicity of the results 
of W7 models (Nomoto et al.lll984i) an d input those pa- 
rame ters into their PHOENIX code (|Hauschildt et al] 
Il996l ) to produce synthetic spectra. The main difference 
between these meth ods is the definition of "metallicity." 
IHoflich et all [l99§) uses the term to mean the metallic- 
ity of the progenitor star, while ILentz et al.l (j20Q0) uses 
it to mean the metallicity of the ejecta. 

The differing definitions of metallicity yield different 
initial conditions, whi ch resulted in contra dictory results 
from these studies. Hoflic h et al.l (|1998f ) suggest that 
with increasing metallicity, the ultraviolet (UV) co ntin- 
uum of the SN increases, while ILentz et~aU ()2000D sug- 
gest that it decreases. Ultimately, the dif ferences are 
the result of differing densit y structures (jLentz et al.1 
[2000'; 'Domingucz et al.' l200fD . Although the method 
of Lcntz e t al. (2000) seems less physical than that of 
IHoflich et al.l (11998) (simply scaling the metallicity of the 
ejecta by solar abundances does not take into account, 
for example, that the Fe-group elements are mainly pro- 
duced in the SN explosion), they provide model spectra 
for varying mctallicitics, which may elucidate differences 
between low and high-redshift SN spectra. 

Further predic tions for lower m etallicity include 
faster ris e times (IHoflich et al.|[l998D . faster light-curve 
decline (IHoflich et al.' '199^, lower ^"^Fe produc- 
tion (Hoflich et al. 1998), sm aller blueshifting of 
Si II A6355 (Lc ntz et al.l |2000D . decreasing B - V 
color (|Domfnguez et al.| 120011 : iPodsiadlowski et aTl 
20061). and chan gi ng lu m inosity (IHoflich et al.1 1 19981 : 
Dommgue z et al.l 120011 : [Podsiadlowski et al.' zOOw; 
Timmes et~aI1 I2003D . |R6pke & Hillcbrandt (2004) 
suggest that the C/0 ratio of the progenitor does not 
signiflcantly affect peak luminosity. 

Observationally, a lack of evolution has been sup- 
ported by investigating various SN quantities such as rise 
time (Riess ct al. 1999a), line velocities (B londin et al.l 
I2006t iGaravini et al.ll2007f) . multi-epoch temporal evo- 
lution (jFolev et al.ll2005f ). line str engths (IGaravini et al.l 
I2007D . and line-strength ratios (|Altavifla et al.l l2006f ). 
There have also been studies comparing the spectra of 
individual hi gh-re d shift SNe la to low-redshift SNe la 
(iRiess et al.l Il9 98: Coil et al. 2000; Ho ok et al.l l2005l : 
iMatheson et al...2005. : .Balland et al.., 20071 ), all of which 
have concluded that there is no clear difference in spec- 
tr al properties between the two samples. 

iBronder et al.l (j2007f ) recently presented measurements 
of line strengths that suggest a difference between low 
and high-redshift SNe la in one of three features mea- 
sured. They find that the difference is highly dependent 
on the galaxy contamination at high redshift and might 
be affected by their small low-redshift SN sample. Con- 
sequently, they note that the difference is interesting but 
not significant. 

Despite the consistencies in spectral properties. 



iHowell et al.l ()2007D note a slight shift in the mean pho- 
tometric properties of SNe la with redshift. They ex- 
plain this evolution as a change in the ratio of pro- 
:enitors from the "prompt" and "delayed" channels 
Scannapiec o & Bildstcn 2005), corresponding to young 
and old progenitor systems at the time of explosion, re- 
spectively. In particula r, the light-c urve shape parame- 
ter "stretch" (fColdhaber et al.ll200i[ ) increases with red- 
shift. HST observations of ESSENCE objects suggested 
that the sample may have a large proportion of objects 
with slow-declining (large stretch) light curves, but this 
is pro bably the result of a selection bias ()Krisciunas et al.l 
[2005h . Since stretch (and other luminosity light-curve 
parameters) is correlated with spectral properties, one 
might expect the spectra of high-redshift SNe la, on av- 
erage, to differ from those of low-redshift SNe la. 

Since all galactic environments at redshift < z < 1.5 
are also present in the local Universe, SN la evolu- 
tion does not necessarily mean that there are not lo- 
cal analogs. For instance, if the distribution of observ- 
ables is on average different at high redshift, as long as 
for each high-redshift SN there is a similar low-redshift 
counterpart, the peak brightness could, in principle, be 
correctly translated into an accurate distance. Within 
the local sample, there is no indication of a correlation 
between host-galaxy metallicity and light-curve shape 
ijCxaflagher et al.ll2005l ). 

In the process of classifying and finding the redshifts 
for SNe from the ESSENCE (Equation of State: SupEr- 
Novae trace Cosmic Expansion) survey (|Miknaitis et al.l 
[2007; Wood-Vasey ct al. 200^, we have obtained 107 
spectra which have accurate light-curve parameters such 
as A (a light-curve width para meter), time of maxi- 
mum light, and vi sual extinction (jMatheson et al.1 [20051 : 
iFolev et al.ll2008bl ) . Most spectra in this sample have low 
signal-to-noise ratios (S/N) compared to spectra of low- 
redshift SNe. This makes impractical a detailed analysis 
of each object individually to test for outliers. However, 
by combining the spectra to make composite spectra, we 
are able to study the mean spectral properties of the 
samples. 

In Section[2|we discuss our low and high-redshift SN la 
spectral samples. We describe our methods of creat- 
ing composite spectra in Section [3] In Section [4| we 
present the composite spectra and compare the two sam- 
ples, while in Section [5| we discuss the implications of 
these results. We present our conclusions in Section [6l 
Throughout this paper we assume the standard cosmo- 
logical model with {h,nm,^A) = (0.7,0.3,0.7). 

2. SUPERNOVA SAMPLES 

In order to test for potential evolution in SN la spectra, 
we need to explo re the largest redshift range possible. 
iRiess et all (|2007() have obtained spectra of 10 objects 
with z > 1, with one SN la at z = 1.39. The large 
look-back time of these objects allows significant time for 
progenitors to evolve between z ~ 1 and 2 = 0; however, 
even at z > 1, many galaxies had already become metal 
rich (e.g.. lCimatti et al.l[200l . 

In addition to these high-redshift SN spectra, through 
the first four years of the ESSENCE campaign we have 
obtained 107 spectra of 88 SNe la with lig ht curves that 
could be fit by MLCS2k2 (|Jha et al.ll2007f) . The spectra 
were obtained by the Keck I and II 10 m telescopes, the 
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Very Large Telescope (VLT) 8 m, the Gemini North and 
South 8 m telescopes, the Magellan Baade and Clay 6.5 m 
telescopes, the MMT 6.5 m telescope, and the Tilling- 
hast 1.5 m telescope at the F. L. Whip ple Observatory. 
The spectr a are analyz ed individually by Mat heson et al.l 
(2005) and lFolev etaD (l2008b'l. whi le the light curves are 
presented bv lMiknaitis et all (|2007f ). The spectra are in 
the redshift range of 0.155 < z < 0.777, with correspond- 
ing rest-frame phases of —10.8 < t < 20.9 d relative to 
i?-band maximum. The objects span a wide range of 
light-curve width from —0.601 < A < 0.843, correspond- 
ing to luminosities of —19.7 < My < —18.6 mag at max- 
imum brightness, respectively. 

The observed wavelength ranges for the objects vary 
because the spectra were obtained with different instru- 
ments and each spectrum has been trimmed individually 
to remove the noisy ends of the spectrum. The rest- 
wavelength ranges are even more disparate because of 
the significantly different redshifts. The rest-frame spec- 
tra span the wavelength range 1940-8174 A. 

For the purpose of comparison, we have also com- 
posed a sample of 397 spectra obtained mainly with 
the Kast double spectrogr aph (iMil ler fc S tond [l993h 
on the Lick 3-m telescope (jFolev et al. 2008^, and 45 
from HST/IUE (to sample the UV portion of the spec- 
trum; SNe 19 80N, 1981B, 198 2B, 1983G, 1986G , 1989B 
Panagiall2007l 1990N. 1991T.ljefev et al.l ll99a 1992A 
Kirshner et al] Il993l 2001eh, and 2001ep ISauer et all 
20081 ). Th e light curves of t hese objects are presented 
elsewhere fButa et al.' '1985: Younger & van den Berghl 
1985: Ciatti et al. 1988; Hamuv et al. 1996; Ricss et ay 
1999bl: IJha etralir2006l : iGaneshahngam et al.ll2008f) . The 
Lick objects are observed an d reduced in a rnanner simi- 
lar to the ESSENC E objects (jMatheson et al.ll2001l . l2005l : 
iFolcv ct al.l l2008bf) . These similarities reduce the system- 
atic differences between observing programs. Although 
the HST/IUE spectra are not necessarily free of reduc- 
tion errors, the differences should be minor and these 
spectra have great utility in our study. 

In Figures [1] [21 and [Sj we present histograms show- 
ing the A, phase, and redshift distributions, respectively. 
The ESSENCE and Lick samples have similar distribu- 
tions, and more importantly, the Lick spectra sample the 
entire range of phases and almost the entire A range of 
the ESSENCE sample, allowing for an unbiased compar- 
ison. The total numbers of spectra in various redshift, 
phase, and A bins are presented in Tables [1] and [2l All 
composite spectra presented in this paper are publicly 
available^ ^. 

3. METHOD 

3.1. Reprocessing 

To properly combine our sample of SN spectra, we 
must first attempt to remove any effects of non-SN 
sources from each spectrum. The two main effects are 
reddening by host-galaxy dust and galaxy-light contam- 
ination. 

The reddening of the spectra causes the spectral shape 
and relative feature strengths to change. To correct for 
the distortion introduced by reddening, we have used the 
MLCS reddening para meters listed in Tabl eland a stan- 
dard reddening law (jCardelh et al.l I1989D to deredden 

|http: //astro .berkeley . edu/$\slm$rf oley/composlte/| 
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Fig. 1. — A histogram of the number of spectra per A bin. The 
thick line with shading represents the ESSENCE objects, while the 
thin line represents the Lick objects. The dotted lines represent 
the distinction between the luminosity subclasses of SNe la (over- 
luminous, normal, and underluminous correspond to A < —0.15, 
—0.15 < A < 0.3, and A > 0.3, respectively) from IJha et all 
igOOi). 
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Fig. 2. — A histogram of the number of spectra per phase bin 
relative to the B-band maximum bin. The thick line with shading 
represents the ESSENCE objects, while the thin line represents the 
Lick objects. 
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Fig. 3. — A histogram of the number of ESSENCE spectra per 
redshift bin. 

our spectra. Since the ESSENCE photometry consists 
of only two bands, it is not possible to independently 
determine both Ay and Ry- To that extent, our cor- 
rected spectra may not be properly dereddened. The rel- 
ative flux difference between a spectrum with Ay — 0.3 
mag with Ry = 2.1 and Ry — 3.1 is less than 10% 
for wavelengths longer than 3400 A, and it is possible 
that the reddening parameters account fo r some intrin- 
sic co lor variability within the SN sample (jConlev et al.l 
l2007f ). However, we have performed our analysis with 
both a low-extinction sample [Ay < 0.3 mag ) and our 
complete, treated sample (see Section yielding sim- 
ilar results. Therefore, the errors in Ay and Ry are not 
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TABLE 1 

Number of Spectra For Redshift/Phase Binning 



Redshift ll<t<-3d -3<i<3d 3<t<10d 10<t<17d 17<t<23d 



53 55 50 46 35 

0-0.2 2 2 

0.2-0.4 7 10 10 8 3 

0.4-0.6 10 8 16 9 2 

0.6-0.8 1 7 8 2 



TABLE 2 

Number of Spectra For Redshift/A Binning 



Redshift 


A < -0.15 


-0.15 < A < 0.3 


A > 0.3 





141 


158 


89 


- 0.2 


2 


2 





0.2-0.4 


15 


19 


4 


0.4-0.6 


26 


21 





0.6-0.8 


15 
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dominant. 

At high redshift, the size of the SN host galaxy and 
the average offset between the host galaxy and the SN 
are small relative to the slits of our spectrographs. Con- 
sequently, it is often difficult to completely remove the 
galaxy-light contribution from the total light at the posi- 
tion of the SN. For some E SSENCE spectra, we employ 
a deconvolution technique (jBlondin et al.l [20051 ). which 
separates the SN and galaxy light better than simple 
background subtraction. With several bands of photom- 
etry for the galaxy, one can properly model the galaxy 
type, allowing one to find a matchin g template spectrum 
to subtract off the SN spectrum (jHowell et alJ l2005f ). 
However, using ESSENCE'S current photometry this is 
not possible. A campaign in underway to image the 
ESSENCE fields in BVz, and the resulting data should 
allow one to properly model the galaxy spectral-energy 
distribution (SED). 

Obtaining a spectrum of the host galaxy after the SN 
has faded and subtracting an appropriate percentage of 
the host-galaxy spectrum from the SN spectrum is the 
best method to remove host-galaxy light contamination. 
This method, however, is typically very time intensive, 
and has not been performed on our sample. We have 
a program underway to obtain the host-galaxy spectra; 
however, observations have not yet begun. In this paper, 
we neglect the galaxy-light contamination in our spectra 
when making the composite spectra and then choose to 
compare our samples in a way that accounts for it (see 
Section HI). 

3.2. Construction of the Composite Spectra 

The composite spectra were constructed by first dered- 
shifting the individual spectra to the rest frame. The 
redshifts have errors up to 0.01 which might artificially 
widen some spectral features in the composite spectrum. 
If we wish to look at a dereddened spectrum, we then 
deredden the spectra based on our values of Ay and Ry 
found from fitting templates to the light curves is 
fixed at 3.1 for all ESSENCE objects, and unless the fits 
require Ry 3.1, the low- re dshift objects also have Ry 
fixed to 3.1; iJha et all 1200711 . We then put the spectra 
through a low-pass filter to remove residual night-sky 



lines, cosmic rays, and host-galaxy absorption and emis- 
sion lines. By comparing the original spectra to a heavily 
smoothed version of each spectrum, we then determine 
the pixel- by-pixel S/N. Finally, we average the spectra, 
weighting by S/N. Since not all spectra have a common 
wavelength range (particularly comparing the lUE spec- 
tra to the Lick spectra), we first constructed a composite 
spectrum from the majority of spectra with overlapping 
wavelengths, and then used that temporary composite 
spectrum (which overlapped with all spectra) to match 
the fluxes of the individual spectra. 

There are some intrinsic issues to constructing a com- 
posite SN spectrum. Despite having hundreds of SN 
spectra, to ensure a sufficient number of spectra per pa- 
rameter bin we must still have somewhat large bins for 
some parameters such as phase. This can smear out cer- 
tain features. However, since SN spectra tend to evolve 
smoothly over small time periods, averaging over a phase 
bin is a reasonable estimate for the average phase spec- 
trum. Other issues are intrinsic differences in SN spec- 
tra. For a given ph ase and A, spectra still differ from one 
object to another llMatheson et al.ll2007T ). Or, within a 
particular bin, objects can have significantly different ex- 
pansion velocities or colors. When constructing a com- 
posite spectrum, most of these spectral differences do 
not change the composite spectrum from looking like the 
"true" average SN spectrum. However, differences in ex- 
pansion velocity will tend to make spectral features wider 
and shallower. Since we are comparing composite spec- 
tra to each other and not a composite spectrum to an 
individual spectrum, as long as the underlying samples 
are similar, the effects of differing expansion velocities 
will not create significant differences in the composite 
spectra. 

The ESSENCE SNe come from a blind search, result- 
ing in many objects in low-luminosity host galaxies. The 
Lick sample, on the other hand, comes primarily from 
targeted searches and is biased to higher-luminosity host 
galaxies. This may result in slightly different samples, 
however the similar A distributions for the ESSENCE 
and Lick samples suggest against this. However, in cre- 
ating composite spectra, since the spectra are weighted 
by S/N, the composite spectra will be weighted more to- 
wards intrinsically overluminous SNe la. Since this effect 
occurs in both samples, this should not affect our analy- 
sis. 

3.3. Determining Spectral Variance 

We want to determine both the average spectrum of 
a given phase. A, and redshift, as well as the variation 
about that average spectrum. To do this, we implement 
a boot-strap sampling (with replac ement) algorithm to 
estimate the variance ()EfrOTllT98l . The variance is a 
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combination of the noise in our spectra, any systematic 
effects during the reduction process (such as poor sky- 
Hne removal and incomplete galaxy-light subtraction), 
and the inherent variance in the SN sample. 

For the Lick sample, the spectra are generally of very 
high quality with little noise or systematic issues. As 
such, the variance in these spectra is dominated by the 
intrinsic scatter among the objects. A detailed analysis 
of these spectra will be presented in a future paper (for 
another low-reds hift composite spe ctrum derived from 
eight objects, see I James et alll2006[ ). For the ESSENCE 
sample, the Poisson noise dominates at the bluest wave- 
lengths since there are fewer spectra adding to the com- 
posite in this wavelength region and the spectra have 
lower S/N in the UV. The reddest wavelengths are dom- 
inated by reduction issues, specifically poor sky subtrac- 
tion residuals. The observed wavelengths correspond- 
ing to night-sky lines are weighted less, but the average 
noise over large wavelength ranges remains higher at ob- 
served near-IR wavelengths than at visual wavelengths. 
The other major factor at long wavelengths is the large 
variation in host-galaxy contamination. Since the SEDs 
will vary dramatically based on the amount of galaxy- 
light contamination (SNe la are blue relative to galax- 
ies), there is additional variance unrelated to the intrinsic 
variability of SN spectra. 

4. RESULTS 
4.1. Comparison of Low-Redshift Templates 

To understand the subtle differences between the 
ESSENCE and low-redshift composite spectra, we must 
first examine the differences between different low- 
redshift composite spectra. 

Th e original Nugent template spectra (iNugent et al.l 
were constructed from 84 spectra (31 for the 
maximum-light template) of which 63% (52% near max- 
imum Hght) come from SNe 1989B, 1992A, and 1994D. 
Two of these objects have been considered slightly atyp- 
ical SNe la , with strong dust absorption (SN 1989B; 
[Wells et al. 1994) or anomalous lurninosity and color s 
(SN 1994D; Richmond et all IT995t IFatat et al.l fl99l . 
Since no SN is the perfect example of a SN la, the pres- 
ence of these objects in the composite is not worrisome. 
However, having the majority of the sample rely on a 
few objects can severely skew the composite away from 
a true average. The Nugent templates have since been 
updated^^, with more objects. However, the influence of 
these atypical SNe la is still strong. 

A new spectral template has been constructed by 
iHsiao et al.l (j2007) for the Supernova Legacy Survey 
(SNLS). This template is created by combining low and 
high-redshift spectroscopy and photometry. It is created 
from a more varied sample than the Nugent template, 
but still has a very limited low-redshift sample. 

The Lick maximum-light composite spectrum, the Nu- 
gent maximum-light template spectrum, and the new 
SNLS template spectrum are shown in Figure [4l The 
Nugent template is bluer at this phase, particularly in 
the UV and the near-UV. If normalized at 4000 A, the 
Lick composite spectrum has higher flux levels than ei- 
ther the SNLS or Nugent template for wavelengths some- 
what redder than 4300 A, although it is very similar to 

http://supernova.lbl.gov/'^nugent/nugent_templates.html . 
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Fig. 4. — {top panel): Composite low-2 maximum-light SN la 
spectra. The black curve is the composite spectrum of the Lick 
sample; the green (dashed) curve is the Nugent template spectrum; 
and the red (dotted) curve is the SNLS template spectrum. The 
grey region is the la boot-strap variation of the Lick composite 
spectrum, (bottom panel): The difference between the Lick com- 
posite spectrum and the Nugent and SNLS template spectra. The 
grey region is the la boot-strap variation of the Lick composite 
spectrum. 

the Nugent template in this region and to both templates 
for wavelengths redder than 4700 A. Besides the feature 
at 4500 A, the SNLS template is very similar to the Lick 
composite. 

The only other region with major differences is in 
the UV (where the Nugent template also deviates from 
the Lick composite spectrum). Both the Nugent and 
SNLS templates are constructed by warping the spectra 
to match the colors of a normal SN la at the phase of 
the spectrum (for the SNLS template, this is done by 
using light-curve information in the construction of the 
template) . Howe ver, the large dispersion in the U band 
pha et al.ll2006[ ) along with extrapolation into the UV 
makes the Nugent template UV continuum dubious. The 
Lick spectra were all observed and reduced in a consistent 
manner, pr oducing relative spectr ophotometry accurate 
to ~±5% ()Matheson et al.ll2000D . The colors for the 
three spectra as well as the colors of the MLCS A = 
template light curves are presented in Table [3l 

Considering the relative diversity of the Lick and Nu- 
gent samples and the methods of producing the compos- 
ite spectra, we consider the Lick composite spectrum to 
be more reliable than the Nugent template. The SNLS 
template uses high-redshift SN data (both light curves 
and spectra) in the construction of their template. As a 
result, this template should not be used for comparison 
to high-redshift SNe la. In particular, the UV portion of 
the SNLS template is heavily weighted to high-redshift 
SNe la. Most importantly, the spectra in the Lick and 
ESSENCE samples were reduced and the composite spec- 
tra were constructed in the same manner. For our pur- 
poses, the Lick composite spectrum is superior to the 
alternatives and will be used for comparison in the rest 
of this paper. 

The relative strengths of most absorption lines in the 
Nugent, SNLS, and Lick spectra are similar. The line 
velocities are also similar, with the exception of Si II 
6355, which has a blueshift at maximum light of —12, 300 
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TABLE 3 

Low-Redshift Maximum-Light Colors 



Spectrum 


U - B 


B-V 


V -R 


V -I 


(mag) 


(mag) 


(mag) 


(mag) 


MLCS A = 


-0.47 


-0.07 


0.00 


-0.29 


Lick 


-0.33 


-0.04 


-0.04 


-0.34 


Nugent 


-0.45 


-0.07 


0.00 


-0.35 


SNLS 


-0.41 


-0.06 


-0.02 


-0.26 



km (Nugent template) versus —11, 200 km s~^ (SNLS 
template) and —11,400 km s~^ (Lick composite). An- 
other interesting difference is the weaker feature at 
3000 A in the SNLS template compared to the Nu- 
gent template and Lick composite spectrum. Although 
iHsiao et al.l ()20 07*) show no major difference at this wave- 
length when including or excluding high-redshift spec- 
tra, the final SNLS template contains both high-redshift 
spectra and light-curve information, which may result in 
a bias for the ~3000 A feature (see Section lO)) . 

4.2. Maximum-Light Spectrum 

The nature of the ESSENCE search dictates that 
most objects are observed spectroscopically near max- 
imum light. As seen in Figure [2l most spectra were 
obtained near maximum light (64/107 have rest-frame 
phases within one week of maximum light), with a me- 
dian of 4.1 d after maximum light. Since objects are 
brightest at maximum light, the highest S/N spectra are 
usually obtained at this phase and the contrast between 
the SN and the underlying galaxy is largest. These fac- 
tors together make the maximum-light composite spec- 
trum a higher S/N spectrum. 

In Figure O we present the ESSENCE maximum-light 
composite spectrum with average parameters (z) — 0.37, 
{t) — —0.8 d, and (A) = 0.01. Comparing the spectrum 
to the Lick maximum-light spectrum (which has average 
parameters (z) = 0.02, {t) = 0.5 d, and (A) = -0.02), 
it is obvious that the two spectra have rather different 
colors, with the ESSENCE spectrum being redder than 
the Lick spectrum. 

As described in Section [44l the large projected size of 
the slit for high-redshift objects makes isolating the SN 
from the host galaxy difficult. It is much easier to sep- 
arate the SN from the galaxy light at low redshift. Our 
ability to remove host-galaxy light from any given SN 
spectrum hinges on how well we can model the galaxy 
background, which is highly dependent on how isolated 
the SN is from its host galaxy. To account for the dif- 
ference in the amount of galaxy light remaining in a SN 
spectrum, we fit a combination of a galaxy spectrum and 
the Lick composite spectrum to m atch the ESSENCE 
spectrum. Similar to the method of (|Howell et al.ll2005l ). 
using five galaxy spectra templates (E, SO, Sa, Sb, and 
Sc) and varying the galaxy light from 0% to 100% of the 
comparison spectrum, we find a best-fit combination. As 
seen in FigureEJ by adding some galaxy light (35% of an 
Sb galaxy spectrum for the maximum-light spectrum) to 
the Lick spectrum, we can match the Lick and ESSENCE 
spectra quite well. We believe galaxy-light contamina- 
tion to be the main factor in the discrepancy between 
the Lick and ESSENCE continua. 

Properly determining the galaxy contamination is cru- 
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Fig. 5. — (top panel): ESSENCE maximum-light composite 
SN la spectrum. The black, relatively noisy curve is the ESSENCE 
composite spectrum. The purple curve, slightly below the black 
curve, is the Lick maximum-light composite spectrum. The red 
(dashed) curve is an Sb galaxy spectrum. The green curve is the 
addition of the red and purple curves, fit to match the ESSENCE 
composite spectrum. The combination is 35% of host-galaxy spec- 
trum in the region 4000—6000 A. The grey region is the Icr boot- 
strap variation of the ESSENCE composite spectrum. The large 
variance at red wavelengths is due to strong sky-subtraction resid- 
uals in individual spectra and not intrinsic variability in the SN 
sample, (bottom panel): The difference between the ESSENCE 
and Lick composite spectra. The grey region is the Itr boot-strap 
variation of the ESSENCE composite spectrum. 

cial in determining differences between our two samples. 
We have performed several tests to determine the validity 
of our claims with different galaxy contamination. 

In addition to fitting template galaxy spectra to our ob- 
served spectra, we have reconstructed the best-fit galaxy 
spectra from galaxy eigenspectra. Using the first four 
Sloan Digital Sk y Survey (SDSS) galaxy eigenspectra 
(|YiD et al.l[200l . we fit the eigenspectra to the resid- 
uals of the ESSENCE and Lick spectra. Higher-order 
eigenspectra are dominated by emission lines and high- 
frequency modulations. 

The process consists of reducing the of 

/essence - (a/Lick + bUgj3gg) , (1) 

where /essence is the fiux vector of the ESSENCE com- 
posite spectrum, /Lick is the fiux vector of the Lick com- 
posite spectrum, Usdss = 1 < i < 4 is the matrix 
of SDSS eigenspectra, a is a free parameter, and b is 
the eigenvector that best describes the /essence " a/Lick 
residual. The best-fit reconstructed galaxy spectrum is 
comprised of 83.5%, 7.2%, 5.6%, and 3.7% of the first- 
four eigenspectra, respectively. The first eigenspectrum 
is the average of the SDSS galaxy sample and resembles 
an Sb galaxy. 

In Figure [6l we present the principal-components anal- 
ysis (PCA) reconstructed best fit, and our Sb galaxy 
template spectrum. The PCA-reconstructed and tem- 
plate spectra are fit with differing amounts of SN light, 
so their absolute scaling is approximately the same, but 
slightly different. Despite this, we can still see that the 
PCA and Sb template spectra have very similar SEDs, 
differing the most around 4200 A and at wavelengths 
longer than 6500 A. Extrapolating the PCA spectrum, it 
appears that the spectra do not differ drastically in the 
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Fig. 6. — The best-fit PCA-reconstructed and template galaxy 
spectra (black and red curves, respectively). Since the amount of 
galaxy light in the fit changes based on the SED of the galaxy 
spectrum, the scaling between the spectra is arbitrary. 



uv. 

Figure [7] shows the fit and residuals of the PCA- 
reconstructed galaxy spectrum (similar to Figure^ ■ The 
residuals for the PCA-reconstructed fit are smaller than 
those for the template fit. Since the addition of the PCA- 
reconstructed galaxy spectrum to the Lick composite 
spectrum gives, by construction, the smallest deviation 
from the ESSENCE composite spectrum, we consider 
any additional differences to come from sources other 
than galaxy contamination. Considering the overall sim- 
ilar results and smaller wavelength range of the PCA- 
reconstructed galaxy spectrum, we will use the template 
galaxy spectrum for comparisons, but note differences 
with the PCA-reconstructed galaxy spectrum when ap- 
propriate. 

It is unlikely that miscalculated reddening corrections 
can cause the differences in the continua. Using the 
unreddened sample of spectra produces a similar contin- 
uum (see Section 14. 3p . Although it is possible that SN 
evolution could cause the discrepancy, no current models 
predict that the spectrum would change in this way. 

After correcting for galaxy-light contamination, by 
adding galaxy light to the Lick composite spectrum^^, 
the Lick spectrum is very similar to the ESSENCE com- 
posite spectrum. Despite the similarities, there are two 
differences worth noting: the ESSENCE composite spec- 
trum lacks the absorption at ~3000 A and has a weaker 
feature at ^^4900 A. Th e feature at ^3000 A is attribute d 
to multiple Fe II lines ([Branch fc Venkatakrishnalll986D . 
while the feature at ~4900 A is attributed to Si II, Fe II, 
and Fe III, with the red wing, where the discrepancy oc- 
curs, being Fe I II A51 29. We will use the convention of 
iGaravini et al.l (|2007f ) and call the entire feature "Fe II 
4800," naming the lines individually only when we are 
discussing particular species contributing to t he fea ture. 

Using the prescription of IGaravini et al.l (|2007f ). we 
measure the pseudo-equivalent widths (pEWs) of several 

Throughout this paper, in all comparisons, ESSENCE com- 
posite spectra are presented unchanged, and Lick composite spec- 
tra have additional galaxy light to match the ESSENCE spectra, 
unless specifically noted. 
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Fig. 7. — (top panel): ESSENCE maximum-light compos- 
ite SN la spectrum. The black, relatively noisy curve is the 
ESSENCE composite spectrum. The purple curve, slightly be- 
low the black curve, is the Lick maximum-light composite spec- 
trum. The red (dashed) curve is the sigma-clipped SDSS PCA- 
reconstructed galaxy spectrum. The green curve is the addition of 
the red and purple curves, fit to match the ESSENCE composite 
spectrum. The combination is 40% of galaxy spectrum in the re- 
gion 4000-6000 A. The grey region is the Icr boot-strap variation 
of the ESSENCE composite spectrum. The large variance at red 
wavelengths is due to strong sky-subtraction residuals in individual 
spectra and not to intrinsic variability in the SN sample, {bottom 
panel): The difference between the ESSENCE and Lick compos- 
ite spectra. The grey region is the Icr boot-strap variation of the 
ESSENCE composite spectrum. 



features in both the ESSENCE and Lick composite spec- 
tra. These values are presented in Table HI We note that 
because of the host-galaxy contamination, the pEWs of 
the ESSENCE composite spectrum are lower limits. We 
have calculated the pEWs of the Lick composite spec- 
trum both with and without galaxy light added. The 
values of the spectrum with galaxy light added should 
be compared directly to the values of the ESSENCE 
spectrum; however, the values without galaxy light give 
both an accurate measurement and an upper limit for 
the pEWs. There is significant systematic uncertainty 
to these measurements from the placement of the con- 
tinuum. We did not attempt to model th e exact system- 
atic e rrors, but they are typically ~10 A ([Garavini et al.l 
[200l . 

For most features, the ES SENCE and Lick com posite 
spectra have similar pEWs. iBronder et al.l (|2007f ) found 
a possible difference in the strength of the Mg II 4300 fea- 
ture, but we see no difference in our com posite spectra. 
We will show a full analy sis similar to the Bronder et aD 
(|2007) study in lFolev et al . (2008b). The Si II 6150 pEW 
has a slightly smaller value in the ESSENCE spectrum. 
Similarly, the Fe II 4800 feature is much weaker in the 
ESSENCE spectrum, which is likely because of a weaker 
Fe III A5129 line. Using the PCA-reconstructed galaxy 
spectrum, we see that the Lick Si II feature has a sim- 
ilar pEW to the ESSENCE composite spectrum. The 
Fe II 4800 feature continues to have a larger (although 
slightly smaller than using the galaxy template spectrum, 
pEW = 94.9 ± 0.2 A) pEW than the ESSENCE compos- 
ite spectrum. The PCA-reconstructed galaxy spectrum 
does not cover the 3000 A region. 
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TABLE 4 
Pseudo-Equivalent Widths 



Spectrum 


Ca II H&K 


Si II 4000 


Mg II 4300 


Fe II 4800 


SII w 


Si II 6150 




(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


Lick 


112.8 (0.6) 


12.4 (0.6) 


86.4 (0.4) 


140.9 (0.4) 


68.5 (0.4) 


99.5 (0.6) 


Lick + Galaxy 


100.7 (0.5) 


9.6 (0.4) 


73.6 (0.3) 


98.8 (0.2) 


43.4 (0.2) 


61.4 (0.3) 


ESSENCE 


99.7 (1.1) 


9.7 (0.5) 


73.8 (0.8) 


78.5 (0.8) 


45.7 (0.8) 


51.7 (1.7) 
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Fig. 8. — ESSENCE low-Ay maximum-light composite SN la 
spectrum. The black curve is the ESSENCE low-Ay compos- 
ite spectrum and the green curve is the ESSENCE full-sample 
maximum-light composite spectrum (as shown in Figure O. The 
spectra contributing to the full-sample spectrum have been dered- 
dened by the values listed in Table [5] The low- Ay spectrum has 
not been dereddened. The low-Ay and full-sample spectra are 
nearly identical, indicating that the difference between the Lick 
and ESSENCE spectra is not the result of miscalculated reddening 
corrections. 

4.3. Low- Extinction Spectrum 

As mentioned in Section 14.21 a possible explanation 
for the difference between the Lick and ESSENCE com- 
posite spectra is a miscalculation of the reddening from 
the light curves. To test this hypothesis, we have created 
an un-dereddened composite spectrum using only objects 
with Ay < 0.3 mag. This subsample is less dependent 
on the precise values of Ay and Rv, reducing those po- 
tential sources of systematic error. Figure [8] shows both 
the full, dereddened and low-Ay (with no dereddening) 
ESSENCE composite spectra. Both the low- Ay and full, 
dereddened ESSENCE samples produce nearly identical 
spectra. There is a slight difference between the con- 
tinua, with the low-Ay sample being slightly redder. 
This is probably due to the slight extinction that has 
not been corrected in the contributing spectra. 

4.4. Redshift- Binned Spectra 

We have divided the full sample of ESSENCE spec- 
tra into four redshift bins with Az = 0.2. Since the 
ESSENCE sample is confi ned to z < 1, we also utilize 
the composite spectrum of iRiess et al.l ()2007D , which has 
(z) ~ 1.1. The composite spectra with mean phases 
(i) = -2.2, -0.6, -0.1, and -0.1 d for the < z < 0.2, 
Q.2 < z < 0.4, 0.4 < z < 0.6, and 0.6 < z < 0.8 bins 
(respectively) are presented in Figure O Broadly, the 
spectra resemble each other with no major differences 



except for redder continua with increasing redshift. We 
attribute the redder colors with increasing redshift to an 
increase in the galaxy fraction of the spectrum. 

Using the method described in Section [42l we deter- 
mined the galaxy-light fraction of each redshift-binned 
composite spectrum. Figure [10] shows the galaxy fraction 
of the redshift-binned ESSENCE composite spectra. For 
comparison, we also show the relationship between the 
projected distance of 1" and redshift, given by 



X dj 



(1 + zf kpc. 



(2) 



where dL is the luminosity distance given in Gpc. If the 
separation between the SN and the galaxy is less than the 
larger of the slit width and the seeing, both typically 1", 
there will be galaxy contamination in the SN spectrum. 
If the galaxy is larger than the seeing, the separation can 
be larger and the SN spectrum can still be contaminated 
by galaxy light. If the physical distance between SN 
and galaxy nucleus, galaxy luminosity, and galaxy size 
do not change, either by galaxy evolution or selection ef- 
fects, we expect the galaxy fraction to increase with d^ 
since larger redshift corresponds to more of the galaxy in 
the slit. Figure [To] shows this trend, further suggesting 
that the main difference between the Lick and ESSENCE 
composite spectra is the result of galaxy-light contami- 
nation. The HST composite spectrum is not affected by 
galaxy-light contamination since the high angular resolu- 
tion and narrow slits of HST give comparable conditions 
to those of the low-redshift objects. Confirming this hy- 
pothesis, our galaxy-light fitting routine yields a best fit 
of 0% galaxy light for the HST composite spectrum. 

In the three lowest- redshift bins, we are able to exam- 
ine the Si II A6355 line. In these bins, it does not appear 
to be stronger at higher redshift, consistent with that 
found for the maximum- light spectrum (see Section 14.21 
and Table H]). The major differences between the to- 
tal maximum-light composite ESSENCE spectrum and 
the Lick composite spectrum are the lack of absorption 
at -3000 A in the ESSENCE spectrum and the weaker 
Fe III A5129 line. Both lines are shown in detail in Fig- 
ure [9] The ^^3000 A line is only available in the three 
highest-redshift bins. In the 0.2 < z < 0.4 bin, there 
appears to be some absorption at this wavelength. The 
0.4 < z < 0.6 bin does not appear to have any absorp- 
tion, but the line does have a positive pEW that is 3.3cr 
from zero. For the highest-redshift bin (z > 0.6), the 
absorption is lacking. The Fe III A5129 line is within the 
wavelength range of all ESSENCE redshift-binned com- 
posite spectra. The lowest two redshift bins appear to 
have line strengths consistent with the Lick composite 
spectrum. However, the highest two redshift bins appear 
to be lacking some absorption. 

The UV portion of the Lick composite spectrum con- 
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Fig. 9. — ESSENCE maximum-light composite SN la spectra 
for different redshift bins. The spectra have average redshifts of 
0.16, 0.33, 0.48, and 0.64, as well as average phases of —2.2 d, 
—0.6 d, —0.1 d, and —0.1 d, respectively. The grey regions are the 
1(T boot-strap variation. The green curves are the Lick maximum- 
light composite spectrum with 34%, 35%, 41%, and 49% galaxy 
light for the comparisons to the (z) = 0.16, (z) = 0.33, (z) = 0.48, 
and (z) = 0.65 b i ns, re spectively. The HST composite spectrum 
from lRiess et all II2007I) . with (z) = 1.1, is composed of individ- 
ual spectra spanning phases of 0—10 d past maximum brightness. 
For comparison, a Lick composite spectrum was constructed from 
low-redshift spectra with similar phases. For the HST composite 
spectrum, no galaxy light was added to the Lick comparison spec- 
trum. We argue that galaxy-light contamination should be lower 
in the HST spectra (see text), and this is confirmed by our galaxy- 
light fitting routine which gives a best fit with no galaxy light. The 
left and right panels show the regions surrounding the ~3000 A and 
Fe II 4800 features. 




Redshift 

Fig. 10. — The galaxy fraction for ESSENCE maximum- 
light composite spectra of varying redshifts found by fitting the 
Lick maximum-light composite spectrum and a varying amount of 
galaxy light to the ESSENCE spectra. The solid curve represents 
the projected size of 1" with redshift. The typical size of a slit 
used to obtain the spectra is 1". 



sists of only a few spectra from even fewer objects. As 
such, the ^3000 A feature might be anomalous for typ- 
ical SNe la. However, the presence of this line in one of 
our ESSENCE spectra but not in the two highest-redshift 
bins suggests that the line has some redshift dependence. 
Since we preferentially select objects that have bright UV 
fluxes at higher redshift, this may be a bias in the sam- 
ple. The (z) — 1.1 composite spectrum from lRiess et alj 
(|2007t l does show a shght dip at 3000 A, but the spec- 
trum is rather noisy at A < 3500 A and the dip is narrow 
compared to other features, so it is likely insignificant. 
If the line is present in the HST spectrum but not in 
the highest-redshift bins of the ESSENCE sample, then 
it is likely associated with another property that changes 
with redshift in the ESSENCE sample, but is unaffected 
in the HST objects, presumably the result of a selection 
effect. Unfortunately, the HST composite spectrum does 
not contain the Fe III A5129 feature, so we are unable to 
confirm the trend to higher redshifts and to determine if 
the trend is the result of a selection effect. 

4.5. Phase-Binned Spectra 

Another subsampling we created was based on phase, 
as estimated from the light curves. We present five 
phase-binned composite spectra in Figure lll[ corre- 
sponding to roughly one week before maximum light, 
maximum light, and one, two, and three weeks after 
maximum light, having phase bins of— ll<i<— 3d 
{{t) = -6.2 d), -3 < t < 3 d {{t) = -0.8 d), 3 < f < 10 d 
{{t) = 5.6 d), 10 < t < 17 d {{t) = 13.3 d), and 
17 < t < 23 d {{t) = 17.9 d), respectively. Again, aU 
ESSENCE composites are generally similar to the Lick 
composites with some differences. 

The -6.2 d ESSENCE spectrum displays three note- 
worthy differences: a higher Si II A6355 velocity (12,800 
km s""'^ versus 12,100 km s~^ for the Lick spectrum), a 
slightly stronger Si II A6355 line, and a UV excess. The 
ESSENCE spectrum has a slightly older average phase 
(—6.2 d compared to —6.6 d) and slightly smaller av- 
erage A (—0.14 compared to —0.08). Although these 
differences are small, they may account for part of the 
different Si velocity and line strength. For some SNe, our 
redshift errors can be as large as ^3000 km s~^ , which is 
much larger than the difference between the two velocity 
measurements. However, the other lines in the spectrum, 
including the Si II A4130 line, do not show this velocity 
shift, indicating that this is intrinsic to the Si II A6355 
line. The UV excess is very difficult to explain as merely 
a systematic error, although see Section [5.21 for a discus- 
sion of the possible issues regarding the treatment of the 
galaxy-light contamination. 

The —0.8 d maximum-light ESSENCE spectrum is dis- 
cussed in Section [42l 

The 5.6 d ESSENCE spectrum is similar to the compa- 
rable Lick spectrum, but there are more differences than 
in the earlier spectra. The blue wing of the S II doublet 
(A5612) is stronger and the overall continuum match- 
ing is worse than with other spectra. This is likely the 
result of poor galaxy-contamination correction. Later- 
epoch spectra require starburst-galaxy spectra for the 
Lick spectra to match the continua of the ESSENCE 
spectra. The 5.6 d spectrum may be the transition be- 
tween these two regimes. 
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Fig. 11. — ESSENCE composite SN la spectra for different phase 
bins. The spectra have phase bins of —11 < t < 3 d, —3 < t < 3 d, 
3 < t < 10 d, 10 < i < 17 d, and 17 < t < 23 d, with average 
phases of —6.2 d, —0.7 d, 5.7 d, 13.3 d, and 17.9 d, respectively. 
The — 3 < t < 3 d spectrum is the same as in Figure |5] The grey 
regions are the la boot-strap variation. The green hues are Lick 
composite spectra with the same phase bins and similar average 
phases. 

The 13.3 d and 17.9 d ESSENCE spectra are also 
similar to their Lick counterparts. However, to make 
this match, it was necessary to include a percentage of 
starburst-galaxy spectra. Restricting our fitting method 
to normal galaxies, all features match between the 
ESSENCE and Lick spectra, but their strengths are very 
different, with the Lick features being much stronger. 
With only normal galaxy light, the ESSENCE spectra 
have very large blue excesses. Unfortunately, these bins 
contain lower-redshift SNe, and as such, we are unable 
to probe the UV portion of the spectrum. 

4.6. A-Binned Spectra 

To test the possibility that a certain luminosity-based 
subsample may deviate from the local comparison, we 
se parated the objec ts into three luminosity bins defined 
bv iJha et all (p006h of A < -0.15, -0.15 < A < 0.3, 
and A > 0.3, corresponding to overluminous, normal, 
and underluminous objects, respectively. 

In Figure 1121 we present the luminosity-binned 
maximum-light ESSENCE composite spectra as well as 
the total and luminosity-binned Lick spectra for com- 
parison. The ESSENCE and Lick spectra are all gener- 
ally consistent with both the total and luminosity-binned 
Lick spectra. The total and luminosity-binned spectra 
are similar, with the total and normal- luminosity spectra 
being nearly identical. This indicates that the total com- 
posite spectra of Lick (and to some degree the ESSENCE 
spectra) are a proxy for the normal-luminosity SNe la. 

The normal-luminosity ESSENCE and Lick spectra are 
very similar, with the spectra having minor differences. 
The underluminous spectra are also very similar; the 
main difference is that the Si II A5972 line is slightly 
stronger in the ESSENCE composite. Since the under- 
luminous composite is constructed from only three spec- 
tra, this difference might not persist for a larger sample. 
Although the overluminous spectra are generally simi- 
lar, there are minor differences between the Lick and 
ESSENCE composite spectra. The Fe II 4800 feature. 
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Fig. 12. — ESSENCE maximum-hght composite SN la spectra for 
different A bins. The composite spectra consist of 3 (10), 14 (18), 
and 15 (9) individual spectra with average A of 0.33, (0.43), 0.01 
(—0.05), and —0.32 (—0.28) for the underluminous, normal, and 
overluminous subsamples defined by Jha ct al. (2006) (A > 0.3, 
-0.15 < A < 0.3, and A < -0.15) for the ESSENCE (Lick) 
sample, respectively. All have average redshifts of ~0.3. The grey 
regions are the la boot-strap variation. The green lines are the 
Lick composite comparison spectra. 

which includes the Fe III A5129 line, is weaker in the 
ESSENCE spectrum, practically vanishing. 

5. DISCUSSION 

Section [?] shows that in every instance, the overall 
appearance of the low-redshift Lick and high-redshift 
ESSENCE spectra is very similar, with only a few small 
but notable differences. In multiple cases, the Fe II ab- 
sorption at ~3000 A present in the Lick composite spec- 
tra is missing in the ESSENCE composite spectra. Sim- 
ilarly, the Fe II 4800 feature is weaker in the ESSENCE 
spectra, particularly at higher redshift and smaller A 
(that is, overluminous SNe la). The Si II A6355 line is 
perhaps slightly stronger in the premaximum ESSENCE 
composite spectrum, but weaker in others. In the pre- 
maximum ESSENCE composite spectrum, there is a UV 
excess and the Si II A6355 line is blueshifted more com- 
pared to the Lick composite spectrum. In this section, 
we attempt to explain the physical nature as well as some 
consequences of these differences. 

5.1. Evolution vs. Changing Demographics 

If there is a true spectroscopic difference between low 
and high-redshift SNe la, it is still necessary to deter- 
mine if objects with the same value of A (and hence the 
same light-curve shape, to first order) are changing with 
redshift (evolution) or if the population of the objects is 
changing with redshift (changing demographics). If, for 
instance, SNe with a given value of A are changing with 
redshift as a result of metallicity differences, we would 
call that "evolution." On the other hand, if the metal- 
licity is lower in star-forming galaxies and more SNe la 
occur in star- forming galaxies at high redshift, causing 
the average SN to have lower metallicity, we would call 
that "changing demographics." 

A difference in samples caused by both changing de- 
mographics and evolution are possibilities. Photometric 
information suggests that SN la demographics do change 
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Fig. 13. — 1(T, 2(7, and 3cr confidence contours for the linear fit 
parameters for the ~3000 A feature. The smaller, solid contours 
are for the Lick sample. The larger, dashed contours are for the 
ESSENCE sample. 

with redshift (: Howell et al.ll2007| ): namely there are more 
objects with higher stretch at high redshift. Until we 
compare objects with similar conditions (host galaxies, 
delay times, etc.), we will not be able to distinguish be- 
tween these possibilities. 

5.2. Line Strengths 
5.2.1. Fell 3000 

The Fe II ^3000 A absorption is absent in the normal- 
luminosity ESSENCE composite spectrum, but the fea- 
ture has similar strength in the overluminous ESSENCE 
and Lick composite spectra. Although the Fe II 4800 
feature is also weaker at higher redshift (because of a 
weaker Fe III A5129 line, which is a component of the 
Fe II 4800 feature), the opposite is true, with the line 
being weaker in the overluminous ESSENCE composite 
spectrum but having similar strength in the normal lu- 
minosity ESSENCE and Lick composite spectra. It is 
therefore unlikely that the lack of the Fe II ~3000 A ab- 
sorption is the result of an ionization or density effect 
(although there may be non-LTE effects). Alternatively, 
there may be another species contributing to the feature. 

To determine the significance of the differences in the 
strength of this feature, we measured the line in indi- 
vidual spectra. After adding the appropriate amount 
of galaxy light (measured from the fit to the spectra), 
we measured the pEW of all low-redshift spectra. We 
then fit a line to the pEWs as a function of age for 
— ll<i<lld. We also performed a linear fit for the 
ESSENCE data over the same age range. In Figure [T51 
we show the confidence contours of the fit parameters 
for the Lick and ESSENCE samples. The contours over- 
lap, indicating that the differences in the samples are not 
significant. 

Nevertheless, we do see some indications of differences 
in the strength of the line with redshift. As seen in Fig- 
ure [TU the Lick composite spectrum shows a strong line, 
while the z > 0.6 ESSENCE and HST spectra show 
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Fig. 14. — Psucdo-equivalent widths of the ~3000 A (squares) 
and Fe II 4800 (diamonds) features as a function of cosmic time. 
The highest-redshift bins are consistent with pEW (3000 A) = 0. 
Overplotted are the best-fit Verhulst (1845) curves showing a tran- 
sition from weaker to stronger features at ttr,'MOO = 9-9 Gyr 
{ztr,3000 = 0.34) and ttr,4800 = 9.7 Gyr (24^,4800 = 0.36). The 
pEW errors shown are only statistical; systematic errors are typi- 
cally 10 /AA. The decrease in the Fe II 4800 fea ture may be arti- 
ficially smaller at higher redshift. See Figure [161 for measurements 
of individual spectra. 

no indication of a line. The lower-redshift ESSENCE 
composite spectra also show weaker lines. We note 
that visually, there appears to be a slight absorption in 
the (z) = 0.33 spectrum, but nothing distinct in the 
(z) — 0.48 spectrum. There may be absorption in the 
HST composite spectrum (see Section |4^ and Figure |9]), 
but it is consistent with no line. Considering the system- 
atic uncertainties, it is possible that the data can be fit 
with a straight line and there is no change with redshift. 

In Figure [TH we plot the evolution o f pEW of the_~ 
3000 A feature with cosmic time. Using a I Verhulst ()1845l ) 
function, which is a logistic function originally used to 
describe population growth, we fit the pEWs assuming 
pEW(t = 0) = and pEW{t = 00) = pEW(z = 0), 
resulting in 

^) - 3.54 X 10^ + 3.16 (e0.^3s._i) A, (3) 

which is also shown in Figure 1141 This equation yields 
a transition time, ttj. = 9.9 Gyr, where the feature is at 
half its current strength. This corresponds to Ztr = 0.34. 
However, considering systematic errors, we believe that 
the feature is consistent with having no change over time. 

The analysis of the ~3000 A absorption is hampered 
by a lack of low-redshift UV spectra. Examining indi- 
vidual spectra, we see a range of strength fro m strong 
absorption (SN 1992A; iKirshner et"alTI 19931) to emis- 
sion at this wavelength (SN 1990N; iJefferv et aD [l992). 
All low-redshift objects (including those without light- 
curve information, which are not included in any com- 
posite spectra) show this feature in absorption in at least 
one spectrum. However, because this wavelength can 
show both absorption and emission, a composite spec- 



12 



Foley et al. 



trum from SNe la could produce a lack of absorption 
(similar to the ESSENCE spectra). Without a better 
understanding of this feature and its strength in the low- 
redshift sample, we do not have the ability to make any 
claims of its variation between low and high rcdshift. 

5.2.2. Fein 5129 

The one substantial difference between the low-A 
ESSENCE and Lick composite spectra is the strength of 
the Fe II 4800 feature. The Fe II 4800 feature is a blend 
of many lines, with the strongest being Fe II, Si II, and 
Fe III. Examining Figures [5] and [H we see that the ma- 
jor difference in this feature between the ESSENCE and 
Lick composite spectra is the reddest portion, attributed 
to Fe III A5129. This feature is weaker in both the 
higher-redshift and the low-A binned ESSENCE spectra. 
The average redshift for the low-A composite spectrum 
(z = 0.51) is significantly larger than the total composite 
spectrum (z = 0.37) and the normal-luminosity compos- 
ite spectrum (z = 0.27). Similarly, as seen in Figure [T5l 
the average A decreases with redshift (A — —0.05, 0.20, 
-0.14, and -0.23 for < z < 0.2, 0.2 < z < 0.4, 
0.4 < z < 0.6, and z > 0.6, respectively). 

Note that the Lick and ESSENCE low-A composite 
spectra have the same average A, (A) = —0.27. There- 
fore, the differences are not the effect of comparing differ- 
ent average parameters. However, based on our current 
data set, we cannot determine if (1) on average the Fe III 
A5129 line is weaker in all SNe la with increasing redshift, 
or (2) the feature is weaker in low-A SNe la at high red- 
shift and the increasing influence of low-A objects on the 
sample for increasing redshifts. Because of this ambigu- 
ity, we cannot distinguish between evolution and chang- 
ing demographics. To address this issue, we would need 
a sufficiently large sample of A > —0.15 SNe la at high 
redshift and A < —0.15 objects over the entire redshift 
range. This would allow us to create A — z binned com- 
posite spectra. Unfortunately, our eleven low-A spectra 
and eight A > —0.15, z > 0.4 (two with z > 0.6) spectra 
do not allow this binni ng. 

Fitting the JVcrhulst (jl845D function to the Fe II 4800 
feature, assuming that pEW(i = 0) = pEW(z = 0.65) 
and pEW(t = oo) = pEW(z = 0), we find 

pEW (Fe II 4800) = 4.13 , ^ 5.42 (eO-^^a. _ i) +^^.7 

(4) 

This equation has a transition time of ttr — 9.7 Gyr, 
corresponding to ztr — 0.36. 

As with all equivalent width measurements, one needs 
to examine if the line strength is changing or if the contin- 
uum level is changing. In our case, there are two sources 
of continuum: the strength of the "emission" lines to 
either side of a feature, and the galaxy light. In order 
to match the pEW(Fe II 4800) values for the z > 0.6 
ESSENCE and Lick composite spectra, the galaxy-light 
component of the ESSENCE spectrum would need to be 
underestimated by 44%. It is possible that the "emis- 
sion" features on cither side of the Fe II 4800 feature are 
weaker in the ESSENCE spectra; however, the remark- 
ably similar maximum-light spectra at all points in the 
range 4000-6000 A (except for this feature), as seen in 
Figure El argues against this. 

As we did with the ~3000 A feature, we fit the pEW 
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Fig. 15. — A distribution of ESSENCE objects as a function of 
redshift. The crosses are values for individual ESSENCE objects, 
while the filled diamonds are the values for the redshift-binned 
composite spectra (including the Lick composite spectra). The 
lack of high-A objects at high redshift is the result of selection 
bias. 

measurements from individual spectra as a function of 
time for both samples. The EW measurements for the 
ESSENCE and Lick samples (with galaxy light added 
to best match the ESSENCE sample) are presented in 
Figure [TBI We only look at the data for i < 5 d, which 
is both where the data are linear and where the line is 
dominated by the same absorption features. Linear fits 
were performed for both the Lick and ESSENCE sam- 
ples with the best fits shown in Figure [111 The confi- 
dence contours associated with the linear fits are plotted 
in Figure [T71 The two samples are different at less than 
the la level. We also considered the Lick sample using 
the PCA-reconstructed galaxy spectrum instead of the 
template-galaxy spectrum. As seen in Figure [71 this re- 
sults in a slightly smaller pEW. Regardless of the galaxy 
template used, with any reasonable value of galaxy con- 
tamination the samples are still different at the ~2cr level. 

The galaxy contamination is a difhcult issue. Ideally, 
we would be able to remove the galaxy contamination in 
each spectrum. However, with the low S/N and small 
wavelength range of each individual spectrum and the 
uncertainty in the amount of galaxy contamination, we 
would need to make large assumptions about the over- 
all continuum shape of each SN spectrum. Our method 
of matching the composite spectra avoids these issues 
by comparing high S/N spectra covering a large wave- 
length range. A test of the accuracy of this method is 
to examine the photometric fluxes of the galaxies and 
SNe at maximum light. From our light curves, we are 
able to determine the galaxy contamination in the ob- 
served R band. Then, assuming the SN spectrum is the 
Lick maximum-light composite spectrum, we can deter- 
mine the galaxy contamination in the Fe II 4800 feature. 
Doing this, we find that the galaxy contamination has 
median and mean values of 17-19% and 22-23%, respec- 
tively, regardless of galaxy template. The la upper limit 



Constraining Evolution in SNe la 



13 



250 r 



200 - 



150 - 



100 - 



o 



50 - 



-50 [ 



K^'^-'^Wo +'•" ++ +0^^" 

+ + 



-10 



-5 
Phase (d) 



Fig. 16. — The EW measurements of Fe II 4800 for indiviudal 
spectra in the ESSENCE and Lick spectra. The black crosses are 
for ESSENCE spectra and the green diamonds are for Lick data. 
The size of the points is proportional to their weight. The best- 
fit lines are solid black and green lines for ESSENCE and Lick, 
respectively. The dashed lines are the la errors in the fits. Each 
Lick spectrum was given a weight o{l/N, where A' is the number of 
spectra from that object used in the fit. This weights each object, 
rather than each spectrum, equally. 



100 




Slope 

Fig. 17. — Icr, 2cr, and Sc confidence contours for the linear fit 
parameters for the Fe II 4800 feature. The large, solid contours 
are for the ESSENCE sample and the dashed and dotted contours 
are for the Lick sample. The dashed and dotted contours are from 
the Icr extremes for galaxy contamination using the Sb and Sc 
templates, respectively. The contours from the best-fit galaxy con- 
tamination are not plotted for clarity, but are close to averages 
between the plotted contours. 

is 47-60% for all galaxy templates. Considering that we 
remove some galaxy light from our spectra while reduc- 
ing our data, we regard these values as upper limits. This 
is consistent with what we have found from matching the 
low and high-rcdshift composite spectra. 

A separate way of examining the Fe lines which avoids 
the majority of galaxy contamination is a ratio of the 
depths of the Fe II and Fe III lines. We measured the 
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Fig. 18.— The Fe II/Fe III ratio for spectra in the ESSENCE 
(black crosses) and Lick (green diamonds) samples. 

depths of these lines from a line segment connecting 
the maxima on either side of the feature, similar to the 
method of determin ing the 7^(Si) ratio first presented by 
iNugent et all (|1995[ ). Since the ESSENCE spectra have 
low S/N, finding a true minimum for the Fe III line is 
difficult. We therefore measured the depth of Fe III for 
all spectra at the wavelength of the minimum in the Lick 
maximum-light composite spectrum. If this wavelength 
is slightly off the true minimum, causing the flux of the 
edge of the line to be measured, we would underestimate 
the true Fe III depth, leading to a larger Fe II/Fe III ratio 
for any given spectrum; however, since all spectra were 
treated in the same way, this should not significantly in- 
fluence the results. The values of this ratio are presented 
in Figure [TSl showing that this ratio evolves linearly with 
age, with the ESSENCE sample having a slightly larger 
Fe II/Fe III ratio on average. However, if we perform a 
two-dimensional Kolmogorov-Smirnov test, the samples 
are different at the 95% level, consistent with a ^2a dif- 
ference. 

There are two possible physical explanations for a 
weaker Fe III line in high-redshift SNe la. The first is that 
high-redshift SNe la have lower photosphcric tempera- 
tures (on average) than their low-redshift counterparts. 
The typical photospheric temperature of maximum-light 
SNe la is 10,000 K, which is very close to the tran- 
sition between Fe I I dominance and Fe III dominance 
(|Hatano et al.|[l999l ). A temperature change of 1000 K 
can dramatically change the ratio of the Fe lines. This 
effect has been use d to explain th e strong Fe III lines in 
SN 1991T (Fihppc nkoeraLllTggl ). 

An alternative explanation is tha t the weaker Fe II I 
line is the result of lower metallicity. iHoflich et al.l ()1998l ) 
shows that the physical quantity most sensitive to chang- 
ing metallicity is the ^''Fe production. As metallicity 
decreases, Y^, the electron fractio n, increases. This, in 
turn, causes less ^''Fe production. iHatano et al.l (|2002l ) 
suggests that the distribution of the Fe isotopes within 
the SN ejecta cause Fe II lines to be generated by ^^Fe, 
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while Fe III lines come from ^''Fe. Therefore, the weaker 
Fe III A5129 line may be an indication of l ower metal- 
licity for high-redshift SN la progenitors. ISauer et al.l 
([2008:^ offers a slightly different explanation: as metal- 
licity decreases, there is less backwarming from UV pho- 
tons, which decreases the temperature and significantly 
reduces the Fe III/Fe II ratio. In both models, lower 
metallicity has the same effect: weaker Fe III lines. 

Finally, the red portion of the entire feature may be 
dominated by a species different from Fe III, such as Si II 
or Fe II. Although this is a possibility, it would require 
small blueshifts {v < 5000 km s^^ ) for those lines. 

5.2.3. St II 6355 

A visual inspection of the Si II A6355 line shows that 
the placement of the continuum of the red wing of the 
line is causing the differences between the ESSENCE 
and Lick maximum-light composite spectra. Although 
this may be a real difference, the galaxy light is com- 
parable to the SN component for wavelengths redward 
of 6000 A. Therefore, the comparison of the ESSENCE 
and Lick spectra is highly dependent on the galaxy spec- 
trum. If we examine the feature when using the PCA- 
reconstructed galaxy spectrum to contaminate the Lick 
composite spectrum, the difference is much smaller. Con- 
sidering this, and also that other lines from intermediate- 
mass elements (including Si II A4130) do not show differ- 
ent pEWs between the Lick and ESSENCE spectra, we 
do not believe the Si II difference is significant. We do, 
however, suggest that the differences in the premaximum 
Si II A6355 fine (outlined in Section IO)) are probably real 
and significant. 

5.3. Metallicity 

The two best observables for determining metallicity 
differences in SN la spectra are the velocity and depth of 
the S i II A6355 fine and the UV flux level (|Hoflich et all 
119981: iLentz et al.l l2000t ). As metallicity increases, the 
opacity of the ejecta increases, and consequently the ve- 
locity of the Si II line becomes more blueshifted while 
the line strengthens. As already discussed in Section [U 
the effect of changing metallicity on the UV continuum 
is disputed. Rather than trying to determine the cor- 
rect model, we first note that all models show a change 
in the UV c ontinuum with changing metallicity. In the 
ILentz et al.) (HOOO) model, with increasing metallicity, 
there is more metal-line cooling, creating a redder SED, 
and at the same time, there is more line blanketing from 
metal lines in the blue, absorbing the flux at those wave- 
lengths. The combination of these effects creates a UV 
deficit at high metallicity. Moreover, differing density 
structures can cause these lines to form at smaller radii, 
decreasing the UV flux level wit h increasing metallicity 
(jHoflich et al .111 9981: iDommguez e t al. 200l]). 

As described bv ILentz et al.l ( 200Q) . most metalHcity 
effects are easier to observe at early times. As seen in Sec- 
tion [431 the premaximum ESSENCE spectrum has both 
a stronger Si II A6355 line with a larger velocity, and 
a UV excess, compared to the Lick premaximum com- 
posite spectrum. Unfortunately, the wavelength region 
where we can easily differentiate the models by UV flux 
is at ~2000 A, which is bluew ard of our spectra. How- 
ever, the general trend of the ILentz et al.l (|2000f ) mod- 
els is that the UV flux level decreases with increasing 



metallicity. If this is the case, then the ESSENCE ob- 
jects would have lower metallic ity than the L i ck obj ects. 
However, using the results of iHoflich et al.l (|1998f ) one 
would arrive at the opposite conclusion — namely, the 
ESSENCE objects have hig her metall icity than the Lick 
objects. According to the .Lentz et al.l (2000 ) models, the 
stronger, more-blueshifted Si II A6355 in the ESSENCE 
composite spectrum is suggestive of higher metallicities, 
conflicting with the UV excess. 

The idea of higher metallicity for high-redshift SNe la 
is counterintuitive. However, there is a possible explana- 
tion if we consider the two-channel model for SN la pro- 
genit ors (jMannucci et al.l 120051 : iScannapieco &: BildstenI 
|2005[) . Examining an extreme example, suppose that at 
2 = all SNe la come from a delayed channel and at 
z = 0.5 all SNe la co me from a prompt c hannel (an 
approximation of realitv: FSuUivan et al.ll2006f ). which we 
will assume to have no delay time for simplicity. Since 
the light travel time from z = 0.5 to z = is 5 Gyr, if the 
delay time for the delayed channel is > 5 Gyr, then the 
low-redshift SNe la would come from lower-metallicity 
progenitors. However, the data suggest that the differ- 
ence in delay times is ~2 Gyr (jScannapieco &: BildstenI 
I2OO5I) . 

There are two major possible sources of systematic 
error. First, the shape of the UV continuum is, to 
some degree, dependent on the galaxy spectrum we fit to 
the ESSENCE spectrum. If we have used the incorrect 
galaxy type or an incorrect amount of galaxy light (both 
of which may vary over our wavelength range), our UV 
continuum could be undersubtracted, causing the UV ex- 
cess. However, the UV fiux of galaxies is small compared 
to the optical flux. Hence, if we have undersubtracted 
galaxy light in the UV, we would expect the effect to 
be even stronger in the optical, which it is not. How- 
ever, if our lower-redshift premaximum spectra tended 
to come from SNe in early-type galaxies and our higher- 
redshift premaximum spectra tended to come from SNe 
in late- type galaxies, then the UV portion of the spec- 
trum (which is only visible in the higher-redshift spectra) 
could be contaminated by relatively bluer galaxy light, 
for which we may not correctly account. 

The other possible source of systematic error is the 
scarcity of premaximum low-redshift UV spectra. There 
are 6 spectra contributing to the UV portion of the Lick 
premaximum composite spectrum (1 from SN 1989B, 4 
from SN 1990N, and 1 from SN 1992A). If these objects 
are atypical, having less UV flux than most SNe la at 
early times, the UV "excess" in the ESSENCE spectrum 
would simply be an artifact of comparing to these atyp- 
ical objects. An argument against this is that spectra of 
SNe 1990N and 1992A are included in other Lick com- 
posite spectra with different phase binning (such as the 
maximum-light spectrum), and they do not have signifl- 
cantly depressed UV flux. 

This result requires us to look at the other spectra in 
detail, searching for additional clues of a difference in 
metallicity. As seen in Figures |5l [51 HH and [T^l there 
are no other obvious UV excesses or extraordinary line 
velocities. As discussed in Section 15.2.21 the difference 
in the Fe III line may be the result of differing metallic- 
ity. Since the outer layers of the progenitor are most 
affected by metallicity and these layers are only seen 
at very early times, a noticeable change of metallicity 
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may not be possible for times past maximum brightness. 
iDominguez et al.l ()2001f ) suggests that the B — V color 
at maximum light decreases with increasing metallicity. 
They show that A{B — V) = 0.05 mag for metallicity 
increasing from \og{Z/ZQ) — —10 to —1.7. We mea- 
sure B-V = 0.22 mag and 0.24 mag for the ESSENCE 
and galaxy-contaminated Lick maximum-light composite 
spectra, respectively, suggesting that this effect may be 
present in our spectra. Unfortunately, the galaxy-light 
contamination reduces our ability to measure accurate 
colors from the spectra, so we cannot determine with 
any certainty if this effect is real. 

One place one might expect to see more pronounced 
evolutionary effects is in the low-A objects. Over- 
luminous SNe la are found in star-forming galaxies 
(iHamuv et al.l 120001: lHow ei m200lD. Based on SN rates 
iMannucci et al.l ()2005[ ) and [Scannapieco fc BildstenI 
(120051) have suggested that there exist two channels for 
SNe la: short (~1 Gyr) and long (> 3 Gyr) delays. The 
SN rate from the short-delay channel is proportional to 
star formation. These data suggest that overluminous 
SNe la tend to come from a short-delay progenitor chan- 
nel, from which the progenitor is more biased by galactic 
environment, and thus galactic evolution, than the long- 
delay channel. As star-forming galaxies increase their 
metallicities with time, they imprint that information in 
the SN la progenitors, creating somewhat different spec- 
tra at high and low redshift. 

Comparing the UV spectra of both the ESSENCE and 
Lick low-A composite spectra, we see no UV flux dif- 
ferences in a subsample that should be particularly sen- 
sitive to these changes. However, the low-A sample is 
where the Fe III A5129 line changes the most between low 
and high redshift. Although we cannot say definitively 
if this is the result of evolution or changing sample de- 
mographics, we see that the high-A Lick and ESSENCE 
composite spectra have weak Fe III, with the ESSENCE 
composite spectra weaker than the Lick spectra. If the 
high-A objects come from an old progenitor population, 
the ESSENCE low-A objects from an intermediate-age 
population, and the Lick low-A objects from a young 
population, the Fe III line may be the tracer of the metal- 
licity evolution through these populations. A more ap- 
propriate subsampling, which might yield more definitive 
resul ts, may be bas ed on the host-galaxy star-formation 
rate (jSullivan et al. 2006; Howell et al. 2007.). 

5.4. K-Corrections 

One of the largest systematic errors associated with 
SN distance determinations is the uncertainty in the K- 
corrections. Ideally, one would have a time sequence of 
perfectly flux-calibrated, contamination-free spectra (ei- 
ther for the object or for all template objects) to syn- 
thesize light curves that perfectly match the observed 
bands in the frame of the object. Since getting multiple- 
epoch spectra of numerous low-redshift objects or ob- 
taining multiple-epoch spectra of a single high-redshift 
SN is very time consuming, we implement K-corrections 
to match photometric observations of high-redshift SNe 
to those of low-redshift SNe. If there is a small difference 
between the SED of the high-redshift SN and its low- 
redshift counterpart near the edge of the filter transmis- 
sion function, errors in the magnitudes (and thus overall 
distance) are introduced. 
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Fig. 19. — The difference between maximum-light synthetic rest- 
frame UBVRI and K-corrected synthetic observed-frame ij-band 
and /-band photometry. The black curve is the difference between 
the ESSENCE composite spectrum and the Lick template. The 
green (dashed) and red (dotted) curves are the difference between 
the Lick composite spectrum and the Nugent and SNLS templates, 
respectively. 

One method of determining the K-corrections for a par- 
ticular object consists of using a template low-redshift 
spectrum, warping it to the observed colors of the high- 
redshift SN, and then d etermining the comparable rest - 
frame mag nitudes ( Kim et al.lll996tlNugent et al.ll2002D . 
This method is particularly sensitive to differences be- 
tween SEDs. 

In order to determine the impact of using one partic- 
ular spectral template for K-corrections, we warped the 
Nugent and SNLS templates and the Lick and ESSENCE 
composite spectra to have the same rest-frame colors, 
redshifted the spectra, and extracted synthetic photom- 
etry. The warping is dependent on where the anchor 
points are placed, but this is not critical for our applica- 
tions. We extracted i?-band and /-band magnitudes for 
all three spectra over a redshift range of < 2 < 1. 
Because of its bluer rest-frame wavelength range, the 
ESSENCE spectrum cannot be used for the lowest red- 
shifts. Furthermore, as seen in Figure [5l the large un- 
certainty in the ESSENCE spectrum at red wavelengths 
makes the /-band magnitudes for z < 0.5 unreliable. 

As seen in Figure [121 all four spectra are generally 
consistent to < ±0.02 mag. The main deviation between 
the Lick composite and Nugent template is at z > 0.8, 
where the near-UV is redshifted into the R band. The 
differences between the spectra in the UV (as seen in 
FigureHI create this large difference. The SNLS template 
is more consistent with both the Lick and ESSENCE 
composite spectra. 

This analysis shows that the systematic error associ- 
ated with K-corrections is small but non-negligible, with 
aK ~ 0.01-0.02 mag f or z < 1, consistent with that found 
bvlHsiao et al.l (l2()07f). b ut larger than that assumed by 
IWood-Vasev et al.l (12007'). We also note that the Nugent 
template has large differences compared to both the Lick 
and ESSENCE composite spectra; it should not be used 
for SNe with z > 0.8. 

5.5. Constraining Evolution 
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Fig. 20. — The percent variation in the Lick maximum-light 
composite spectrum. The dashed Hues are the ±10% variation. 
The variation is larger at shorter wavelengths. 




Fig. 21. — The difference in measured i?-band magnitudes for 
the upper and lower la variance Lick maximum-light composite 
spectrum with redshift. This difference indicates an upper limit 
on the i?-band difference after already normalizing the spectra as 
described in Section |3] 



Ultimately, studies of SN evolution are aimed at deter- 
mining how the luminosity of the average SN la changes 
with redshift. The change in luminosity will directly af- 
fect our measurements of SN distances, and thus mea- 
surements of cosmological parameters. Currently, the 
systematic error associated with SN evolution is esti- 
mated to be < 0.02 mag, corresponding to an error of 
0.02 in measuring w, the equation-of-state parameter for 
dark energy. With our composite spectra, we have the 
first way to quantify this error. 

In Figure I20i we show the percent variation (the ra- 
tio of the la error residuals to the SN flux) for the Lick 
maximum-light composite spectrum. For most optical 
wavelengths, these values are dominated by the intrinsic 
variation of the SN spectra. In the UV, however, the 
noise in the composite spectrum from the small num- 
ber of spectra dominates over the variance amongst ob- 
jects. It is noteworthy that for wavelengths in the range 
of 3200-4000 A, the variation is larger than for longer 
wavelengths and increasing with smaller wavelengths. 
This has been seen photometrically as a relatively large 
variation in the ?7-band light curves of nearby objects 
pha et al.ll2006f) . We therefore expect this intrinsic vari- 
ation to continue to increase into the UV. 

We see that the intrinsic variation of SN spectra is at 
~3% for most of the optical range. Binning over larger 
wavelengths, the variation in individiual features will av- 
erage out, causing the extremes of this variation to de- 
crease slightly. In Figure [211 we show the difference in 
observed i?-band magnitudes of the la variance spectra. 
This measurement is an indication of the variation of the 
spectra at maximum light and not the peak luminosity, 
since we have normalized our spectra to have the same 
flux at a particular wavelength. For z > 0.4, we see that 
the difference increases with redshift as the rest-frame 
UV is redshifted into the observed-frame R band. This 
is probably an overestimate of the maximum difference, 
since most spectra will vary within the la range rather 
than being at either the high or low end. 

We also present the percent variation for the 
ESSENCE maximum-light composite spectrum in Fig- 
ure [22l For most of the rest-frame optical, the variation 
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Fig. 22. — The percent variation in the ESSENCE maximum- 
light composite spectrum in grey. The green line is the percent 
variation of the residuals of the Lick maximum-light composite 
spectrum (as seen in Figure [Sjl. The dashed lines are the ±10% 
variation. 

is ^5%. It is possible that the high-redshift sample has 
a larger intrinsic variation than the Lick sample, but it 
is more likely that the variation is dominated by the dif- 
ferent amounts of galaxy contamination and noise in the 
spectra. Figure [221 also shows that the difference between 
the Lick and ESSENCE maximum-light composite spec- 
tra differs by < 10% over most of the rest-frame optical. 
We can therefore constrain evolution between our two 
samples to ~10% across the spectrum. 

6. CONCLUSIONS 

By combining many low-S/N, high-redshift SN la spec- 
tra, we are able to construct the first series of composite 
SN la spectra based on the parameters of redshift, phase, 
A, and Ay- In addition, we constructed similar com- 
posite spectra from a high-quality sample of low-redshift 
SN la spectra obtained over the last two decades. Com- 
parison of the composite spectra has shown that once we 
account for galaxy-light contamination, the two samples 
are remarkably similar. There are several minor devia- 
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tions between low and high-redshift samples. These de- 
viations fall into three categories: related to metallicity, 
related to ^^Ni production, and unknown. 

The UV excess in the ESSENCE prcmaximum spec- 
trum is indicative of a different metallicity for the low and 
high-redshift SNe. Depending on the mode l, a UV excess 
is the result of higher o r lower metallicity (|Hoflich et al.l 
119981 iLentz et al.il2000D . The stronger, more blueshifted 
Si II A6355 line in the ESSEN CE premaximum s pectrum 
indicates a higher metallicity (|Lentz et al.ll2000f) . 

The most significant difference between our samples is 
the varying strength of the Fe III A5129 line. The evo- 
lution of the Fe III line may indicate that high-redshift 
SNe la have lower temperatures than low-redshift SNe la. 
This, in turn, suggests that SNe la should be less lumi- 
nous at high redshift. The weak Fe III line may indicate 
lower ^*Fe production, which could be the result of lower 
metallicity. Alternatively, lower metallicity may cause 
less backwarming from UV photons, decreasing the tem- 
perature and the Fe III/Fe II ratio. It is unclear if this 
difference is from evolution in all SNe la, evolution just 
in the low-A SNe la, changing demographics, or a selec- 
tion effect. Low-A SNe la tend to come from the short- 
delay progenitor channel. Because of the short delay, 
the progenitors of these SNe are more biased by galactic 
environment, and thus galactic evolution, than the pro- 
genitors of long-delay channel SNe la. It is therefore not 
surprising that the difference in the Fe III line is more 
obvious in the low-A objects. 

It is possible that the different strengths of the ^3000 A 
Fe II feature between low and high redshift is an artifact 
of the construction of the composite spectra. An analysis 
of the individual spectra of both samples indicates that 
the samples are not significantly different; however, the 
small number of UV objects hampers this study. 

It is difhcult to definitively detect metallicity differ- 
ences between the Lick and ESSENCE samples. First, 
we have observed three differences between the samples 
which harbinger a difference in metallicity: a UV excess, 
a stronger and more blueshifted Si II A6355 line, and a 
weaker Fe III A5129 line. The UV excess is an ambiguous 
indicator since the models disagree if it indicates lower 
or higher metallicity. The Si II line in the premaximum 
spectrum suggests higher metallicity, and the Fc III line 
suggests lower metallicity. 

We have also shown that the previously published low- 
redshift template spectra have multiple drawbacks when 
comparing to high-redshift composite spectra. We there- 
fore caution against using these templates for studies of 
SN la evolution. Furthermore, deriving K-corrections 
from any low-redshift template is difficult; however, the 
systematic errors are likely to be relatively small. 

We see that the intrinsic variation of low-redshift SN 
spectra is '^3% in the optical. The spectra vary more 
in the near-UV and UV as suggested by photometry 
(|Jha et al.ll2006[ ) . We are able to put the first constraints 
of SN la evolution to < 10%. 

The results of this study are very suggestive, but re- 
quire further investigation. In order to improve our un- 
derstanding of SN la evolution, we propose three future 
studies related to this work. First, the theoretical models 
of the effects of metallicity on SN la spectra should be 
expanded. With the current ambiguity amongst mod- 
els, we cannot determine the direction of the trend in 



metallicity. Second, we should gather many more high- 
redshift spectra to disentangle the redshift-A ambiguity. 
Finally, further UV observations of nearby SNe la are 
desperately needed. The only current instrument avail- 
able for the task is the Swift UVOT. However, previous 
attempts at obtaining SN la U V spectra have been dis- 
appointing (jBrown et al.ll2005l ). We suggest an intense 
campaign spending several hours per spectrum (similar 
to lUE) with the UVOT. In the near future, we may 
once again have the ability to obtain high-quality UV 
spectra with HST using STIS or COS. If the upcoming 
HST servicing mission is successful, we strongly suggest 
a massive campaign to observe local SNe la in the UV. 
Since JWST does not have the capabilities to observe the 
UV, this may be our last opportunity for many years. 
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TABLE 5 
HiGH-2 SN lA Information 



ESSENCE 


lAU 


z 


Phase'' 




Av 


Rest Wavelength 


Name 


Name 




(d) 






Range (A) 


bOlO 


2002iy 


0.591 


-5.0 


-0.166 


0.104 


2517 - 6461 


bOlO 


2002iy 


0.591 


-1.8 


-0.166 


0.104 


2130 - 5832 


bOlO 


2002iy 


0.591 


13.9 


-0.166 


0.104 


3469 - 5933 


bOlO 


2002iy 


0.591 


14.4 


-0.166 


0.104 


2756 - 5605 


b013 


2002iz 


0.426 


-1.0 


0.034 


0.170 


2842 - 7190 


bOlS 


2002iz 


0.426 


19.8 


0.034 


0.170 


3085 - 6255 


b016 


2002ja 


0.329 


0.1 


0.190 


0.359 


2997 - 7725 


b020 


2002jr 


0.425 


-8.4 


0.059 


0.202 


2810 - 7204 


b020 


2002jr 


0.425 


12.6 


0.059 


0.202 


3684 - 6610 


d033 


2003jo 


0.531 


4.8 


-0.322 


0.085 


2779 - 5833 


d058 


2003jj 


0.583 


1.1 


-0.470 


0.119 


2736 - 5641 


d083 


2003jn 


0.333 


3.8 


-0.273 


0.084 


3174 - 6699 


d084 


2003jin 


0.519 


3.6 


-0.197 


0.221 


2821 - 5879 


d085 


2003jv 


0.401 


-0.4 


-0.228 


0.182 


3790 - 6795 


d086 


2003ju 


0.205 


-4.8 


-0.177 


0.628 


4929 - 7983 


d086 


2003ju 


0.205 


20.9 


-0.177 


0.628 


3070 - 8174 


d087 


2003jr 


0.340 


14.7 


-0.112 


0.126 


4044 - 7208 


d089 


2003jl 


0.436 


10.2 


-0.198 


0.134 


2973 - 6218 


d093 


2003js 


0.363 


-2.8 


-0.365 


0.077 


3117 - 6552 


cl42 


2003js 


0.363 


15.0 


-0.365 


0.077 


3111 - 6555 


d097 


2003jt 


0.436 


4.5 


-0.317 


0.116 


2957 - 6218 


d099 


2003ji 


0.211 


18.7 


-0.110 


0.118 


4475 - 7976 


dll7 


2003jw 


0.309 


-4.7 


0.298 


0.209 


3254 - 6822 


dl49 


2003jy 


0.342 


-8.1 


-0.214 


0.170 


3182 - 6654 


e020 


2003kk 


0.159 


-2.4 


-0.034 


0.437 


3192 - 7592 


e029 


2003kl 


0.332 


-1.7 


0.219 


0.259 


2777 - 6606 


el21 


2003kl 


0.332 


0.5 


0.219 


0.259 


3021 - 6906 


el08 


2003km 


0.469 


-10.8 


-0.280 


0.071 


2518 - 6126 


el08 


2003km 


0.469 


-10.1 


-0.280 


0.071 


2566 - 6739 


el08 


2003km 


0.469 


-10.1 


-0.280 


0.071 


2952 - 6072 


el32 


2003kn 


0.239 


-5.7 


-0.128 


0.739 


3410 - 7211 


el36 


2003ko 


0.352 


-1.2 


0.332 


0.304 


3125 - 6609 


el38 


2003kt 


0.612 


6.5 


-0.284 


0.103 


2660 - 5543 


el40 


2003kq 


0.631 


-1.4 


-0.187 


0.121 


2608 - 5478 


el47 


2003kp 


0.645 


1.0 


-0.174 


0.071 


2586 - 5432 


cl48 


2003kr 


0.429 


-7.1 


-0.107 


0.102 


3010 - 6192 


cl49 


2003ks 


0.497 


10.7 


0.016 


0.186 


2873 - 5910 


foil 


20031h 


0.539 


7.9 


-0.090 


0.157 


2067 - 6107 


f041 


20031c 


0.561 


4.8 


-0.301 


0.086 


2144 - 5952 


f076 


200311 


0.410 


3.1 


0.175 


0.227 


2601 - 6259 


f076 


200311 


0.410 


3.1 


0.175 


0.227 


2272 - 6462 


f096 


20031m 


0.412 


3.7 


0.171 


0.324 


2532 - 6575 


1216 


200311 


0.599 


6.2 


-0.104 


0.117 


2512 - 5809 


£231 


20031n 


0.619 


5.7 


-0.247 


0.089 


2193 - 5752 


£235 


20031j 


0.422 


3.4 


0.165 


0.133 


2365 - 6535 


£244 


200311 


0.540 


5.4 


-0.020 


0.131 


2283 - 6064 


£308 




0.394 


0.5 


0.031 


0.145 


2482 - 6271 


g005 


2004£h 


0.218 


2.4 


-0.253 


0.428 


2921 - 7602 


g050 


2003£n 


0.633 


-1.2 


-0.318 


0.128 


3245 - 5835 


g052 


2004£m 


0.383 


-0.5 


0.382 


0.143 


2516 - 6718 


g055 


2004fk 


0.302 


5.1 


-0.294 


1.009 


4086 - 7327 


g097 




0.340 


10.4 


-0.289 


0.322 


2973 - 6791 


gl20 


2004£o 


0.510 


-0.9 


-0.286 


0.186 


2203 - 6160 


gl33 




0.421 


19.5 


-0.351 


0.452 


3743 - 6713 


gl42 




0.399 


13.8 


0.210 


0.523 


2809 - 6719 


gl60 


2004£s 


0.493 


9.8 


-0.308 


0.194 


3576 - 6403 


g240 




0.687 


9.5 


-0.163 


0.062 


1915 - 5498 


h283 


20041ia 


0.502 


1.4 


0.090 


0.265 


2331 - 5872 


h300 




0.687 


8.5 


-0.279 


0.076 


1940 - 5503 


h311 


2004hc 


0.750 


5.9 


-0.491 


0.065 


2018 - 5300 


h319 


2004hd 


0.495 


-5.8 


-0.274 


0.159 


3571 - 6381 


h323 


2004he 


0.603 


-0.9 


-0.108 


0.120 


3331 - 5951 


h342 


2004h£ 


0.421 


11.5 


-0.356 


0.085 


3659 - 7037 


h359 


2004hi 


0.348 


11.1 


-0.182 


0.299 


3030 - 5708 


h363 


2004hh 


0.213 


3.4 


0.O65 


0.775 


4154 - 7666 


h364 


2004hj 


0.344 


5.3 


-0.006 


0.087 


2752 - 7293 


k396 


2004hk 


0.271 


-5.0 


0.843 


0.175 


3202 - 7633 


k425 


2004hl 


0.274 


10.8 


-0.021 


0.250 


2360 - 5767 


k429 


2004hm 


0.181 


-1.3 


-0.094 


0.126 


3234 - 8000 


k430 


20041111 


0.582 


0.7 


-0.112 


0.113 


2261 - 5891 


k441 


20041iq 


0.680 


2.3 


-0.203 


0.091 


3142 - 5678 


k448 


2004hr 


0.401 


0.7 


0.003 


0.311 


2872 - 5888 


k485 


2004hs 


0.416 


5.5 


-0.230 


0.849 


2895 - 6207 


m026 




0.653 


12.0 


-0.201 


0.102 


3361 - 6215 
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TABLE 5 — Continued 



Name 


lAU 

Name 




z 


Phase^ 
(d) 




Av 
(mag) 


Rest Wavclcng 
Range (A) 


m027 




n 




6 


6 


-0 


134 





362 


3383 - 


- 6903 


m027 




n 




9 


1 


-0 


134 





362 


2690 - 


- 7230 


m032 




Q 


155 


15 


4 


-0 


158 





092 


3688 - 


- 8052 


111O34 




n 


ouz 


9 


2 


-0 


087 





130 


2786 - 


- 5684 


m039 




Q 


249 


16 


9 


-0 


223 





507 


3482 - 


- 7106 


m039 




Q 


249 


17 


1 


-0 


223 





507 


4191 - 


- 7833 


m043 




n 


zuu 


11 


4 


-0 


466 


1 


033 


3436 - 


- 7012 


m057 




Q 


184 


14 


1 


-0 


601 





218 


3673 - 


- 7497 


m062 




n 


'^A 7 
OJ- / 


3 


7 





149 





135 


3721 - 


- 7424 


ml38 




Q 


582 


-6 


7 


-0 


306 





073 


2750 - 


- 5612 


ml38 




Q 


582 


-6 


6 


-0 


306 





073 


3128 - 


- 5877 


ml38 




n 


Ooz 


-5 


9 


-0 


306 





073 


2247 - 


- 5897 


ml58 




Q 


463 


10 


7 


-0 


334 





222 


3133 - 


- 4812 


ml58 




n 


^uo 


11 


8 


-0 


334 





222 


3636 - 


- 6520 


ml93 




Q 


341 


6 


9 


-0 


124 





109 


3187 - 


- 5249 


m226 




Q 


671 


13 


9 


-0 


227 





113 


3183 - 


- 5569 


n256 




Q 


631 


6 


1 


-0 


339 





075 


2667 - 


- 5441 


n258 




Q 


522 


6 








032 





110 


2858 - 


- 5831 


n263 




Q 


368 


2 


7 





054 





075 


3181 - 


- 6488 


n278 




Q 


309 


5 


6 





033 





156 


3321 - 


- 6779 


n285 




Q 


528 


13 


5 


-0 


122 





170 


2848 - 


- 5808 


n326 




Q 


268 


-2 


7 





703 





144 


3431 - 


-6999 


n404 







216 


-4 


3 


-0 


069 





694 


3575 - 


- 7298 


p425 







453 


11 


3 





183 





224 


2994 - 


- 6108 


p425 







453 


12 


1 





183 





224 


3423 - 


-6399 


p454 







695 


2 


1 


-0 


276 





070 


2566 - 


- 5236 


p455 







284 


-10 


6 


-0 


061 





229 


3386 - 


- 6911 


p524 







508 


9 


2 


-0 


229 





155 


2824 - 


- 4668 


p528 







777 


4 


1 


-0 


310 





041 


2397 - 


- 3961 


p528 







777 


5 


9 


-0 


310 





041 


2807 - 


- 5228 


p534 







615 


-3 


4 


-0 


096 





107 


3094 - 


- 5980 


p534 







615 


-2 


3 


-0 


096 





107 


3294 - 


- 5736 



Rest-frame phase relative to _B-band maximum.'^ Mv(t = 0) = -19.504 mag-(- 0.736A -|- 0.182A2 + 5 logio(Ho/65) Oha et al.ll2007h . 



